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Abstract 



A combined analysis of the terrestrial neutrino experiments and the 



^1^' Kamiokande observation of atmospheric neutrino anomaly is performed under 



the assumption of the existence of dark- matter-mass neutrinos, as suggested 
by the recent Los Alamos experiment. In the three-flavor mixing scheme of 



X 

H 

I neutrinos it is shown that the constraints from these experiments are so strong 

that the patterns of mass hierarchy and flavor mixing of neutrinos are deter- 
mined almost uniquely depending upon the interpretation of the atmospheric 
neutrino anomaly. 
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There have been accumulating indirect evidences for nonvanishing masses and the flavor 
mixings of neutrinos. They include the solar neutrino deflcit []1| which may be interpreted by 
either the Mikheyev-Smirnov-Wolfenstein mechanism [Q] or the vacuum neutrino oscillation 
0, both being based upon the notion of flavor mixing. The second in the list is the atmo- 
spheric neutrino anomaly flrst observed by the Kamiokande experiment [Q] and subsequently 
conflrmed by other detectors which strongly indicates the large-angle flavor mixing of 
neutrinos. 

The recent announcement of the discovery of nonzero neutrino mass by the Liquid Scin- 
tillator Neutrino Detector (LSND) experiment in Los Alamos may have brought the 
flrst direct evidence for the neutrino masses and the flavor mixing. The experiment may 
have observed the neutrino oscillation — > z7g with oscillation parameters Am^ ^1 — 12 
eV^ and sin^ 29 ~ 10^^ — 10~^'^, if interpreted by the two-flavor mixing scheme. The result 
may be marginally compatible with the earlier results obtained by the Los Alamos and 
the BNL experiments [0] and by the KARMEN collaboration experiment Clearly the 
result, if conflrmed by the continuing runs, has tremendous implications to particle physics 
and cosmology |T2|-|I? . 



In this paper we try to extract the implications of the possible existence of the dark- 
matter- mass neutrinos, as suggested by the LSND result, in the light of the experimental 
informations from the underground, the reactor and the accelerator experiments. We flrst 



observe, as many authors do [^^-|Tj], that one cannot explain the above three phenomena 
simultaneously by the three-flavor mixing scheme without introducing sterile neutrinos. It is 
simply due to the fact that the three-flavor scheme cannot accommodate three hierarchically 
different mass scales, Am^ ~ (1 — 12)eV^ for LSND, Am^ ~ 10~^eV^ for the atmospheric 
neutrino anomaly, and Am^ ~ 10^^ — lO^^eV^ (~ lO^^^eV^) for the MSW (vacuum mixing) 
solution of the solar neutrino problem. 

We derive the constraints imposed on neutrino masses and mixing angles via a combined 
analysis of the reactor and the accelerator data and the atmospheric neutrino anomaly 
under the assumption that at least one of the neutrinos has a mass which falls into the mass 



range 1-10 eV which is appropriate for cosmological hot dark matter. This assumption will 
be referred to as the assumption of dark-matter-mass neutrinos (DMMN) hereafter. We 
employ the mixing scheme based on three-generation neutrinos, as beautifully confirmed by 



the LEP experiments [^]. It will be demonstrated that it is essential to use the three-flavor 
mixing scheme, rather than optional use of various two-flavor mixings, for drawing correct 
interpretation of the data. We will also consider the restrictions imposed by the neutrinoless 
double (3 decay IT^ . 

Amazingly, the constraints imposed by a minimal set of data, the atmospheric and the 
Bugey [T^ experiments, and the assumption of DMMN are so restrictive as to determine the 



masses and the mixing patterns of three flavor neutrinos. Only a few patterns are allowed: 
(A) light 'Vg" and almost degenerate strongly mixed heavy "z/^" and "z/,-", and its mass- 
inverted version, or (B) light 'V^-" and almost degenerate strongly mixed heavy "z/g" and 
'V^", and its mass-inverted one. The choice of the solutions (A) or (B) is dictated by the 
interpretation of the atmospheric neutrino anomaly. The pattern (A) follows if we interpret 
the atmospheric neutrino anomaly as due to the oscillation, while (B) results if it 

is due to the oscillation. 

We stress that neither combinations with atmospheric nor with solar neutrino data are 
compelling. The interpretation is somewhat more involved in the latter case because there 
still exist three types of solutions to the solar neutrino problem based on the neutrino flavor 
mixing; the small and the large-angle MSW solutions in addition to the vacuum mixing 
one. The analysis of this combination will be presented elsewhere ||18[. Nonetheless, we 



should mention that our theoretical prejudice prefers the case with atmospheric neutrino 
data over the other one. It is natural to introduce a sterile neutrinos to accommodate 
the third experimental data left over in both cases. It is, however, difficult to explain the 
atmospheric neutrino anomaly by introducing the mixing with sterile neutrinos. In doing so 



one encounters the trouble with the light-element nucleosynthesis ||19|| . 

We make use of one crucial aspect of the atmospheric neutrino data in our analysis. 
Namely, the Kamiokande group recently provided a new data set called the multi-GeV sam- 



pie [0. They consist of the events with higher energy, > 1.33 GeV, than the previously 
reported data. The important feature of the new data is that, because of the higher en- 
ergy, the path-length dependence of the oscillation probability can be probed by measuring 
the zenith-angle dependence. It is striking that it can be perfectly fitted by the neutrino 
oscillation with mixing parameters Am^ ^ 10^^ eV^ and sin^20 ^ 1 PD[. Such quantita- 



tive agreement with the zenith-angle dependence is the strongest support for the neutrino 
oscillation interpretation of the atmospheric neutrino anomaly. 

We now make three basic observations in view of the data of the LSND and the 
Kamiokande atmospheric neutrino experiments. (1) To have gross deficit in the ratio 
(z/^ + z7^)/(z/e + z7g) we need at least one large mixing angle. (2) To be consistent with 
the rate of the oscillation events of the order of ~ 5 x 10^^ or less at least one mixing angle 
has to be small. (3) The feature of the atmospheric neutrino data indicates that the one of 
the three Am?- is of the order of ~ 10~^eV^. Notice that we cannot have two Am?- of the 



order of 10~ eV because it contradicts with the assumption of DMMN, Am > leV . 

Based on these observations we classify the hierarchy of the neutrino masses into the 
following two types: 

a : mg ^ m2 ^ m^^ b : m^^ ^ mg ~ (1) 

Here the symbols ~ and 3> imply the differences by ~ 10"^eV^ and ~ 1 — lOOeV^, respec- 
tively. Throughout the analysis in this paper the relative magnitude of the masses connected 
by ~ does not matter. The other types of mass hierarchies which are obtained by permuting 
1,2, and 3 will automatically be taken care of because they merely represent relabeling the 
mass eigenstates. 

We recollect the basic formula of the oscillation probabilities with three flavors of neu- 
trinos. We introduce the neutrino mixing matrix U which relates the flavor- and the mass- 
eigenstates as Ua = UaiUi, where the flavor index a runs over e,fi and r and the mass- 
eigenstate index i runs over 1 to 3. We assume the CP invariance in the present analysis. 
In this case the mixing matrix U is real and contains only three angles 6ij. 



With the mixing matrix the oscillation probability of neutrinos of energy E after travers- 
ing the distance L can be written as 



Sap-^i: U^iUpiU^jUp, sm\—^), (2) 

j>i 



(3) 



AE 

where Am?^ = |m? — m^|. As a convenient parametrization of the matrix U we use the 
so called standard form of the Kobayashi-Maskawa matrix, which is now adopted for the 
neutrino mixing matrix: 

C12C13 S12C13 Sl3 

U = -S12C23 - C12S23S13 C12C23 - S12S23-S13 S23C13 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 

where Cij and Sij are the short-hand notations for cos% and sin 9ij, respectively. We note 
that the three real angles can all be made to lie in the first quadrant by an appropriate 
redefinition of neutrino phases. 

The expressions of the oscillation probability are rather cumbersome involving many 
angle factors. Therefore, we shall derive the approximate formulas by taking into account the 
mass hierarchies and the experimental parameters of the three experiments. The oscillation 
probability which corresponds to the LSND experiment is approximately given by 

^(^V l^e) = 4Ci2Ci3(si2C23 + Ci2S23Sl3)^ Sin^(^^L). (4) 

where two terms with Am^2 ~ ^"^13 (which differ only by lO^^eV^) are combined and 
the term with Am23 ignored. The former procedure can be neatly done by utilizing 
the orthogonality relation of the mixing matrix. The latter approximation is completely 
legitimate because the term is smaller than others by factor 10~^ — 10~^ owing to the mass 
hierarchy Am|3/Am^2 — 10~^ ~ 10~^. 

If the atmospheric neutrino anomaly is attributed to the Vj- oscillation the relevant 
formula is 
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P{U^ Vr) = 2(si2C23 + Ci2S23Sl3)^ (Sl2S23 " Ci2C23Si3)^ 

Aui^ L 

+ 4C23S23Ci3(ci2C23 - S12S23S13) (C12S23 + S12C23S13) sin^( )■ (5) 

In (^) the sine-squared factors with large Am^ of > leV^ are replaced by the average value 
i, which can be justified because of the rapid oscillations; the argument of the sine is 
~ 10 - 10^(10^ - 10^) for L= 10(10^) Km for Am^ = 1 - lOOeV^ and E= IGeV 

We note that there exist the possibility that the atmospheric neutrino anomaly is due 
to the z/^ — > z/g oscillation, the possibility one might not naively expect. It is perfectly 
consistent with a small event rate in the LSND experiment because the relevant scales of 
path length and neutrino energy involved in these two experiments are much different. In 
this case the formula for the oscillation probability to be used is 

P{l^fi l^e) = 2Ci2Ci3(si2C23 + Ci2S23Si3)^ 



4Si2S23Ci3Si3(ci2C23 - S12S23S13) sin^( 777—)- (6) 



Finally the formula for the Bugey experiment takes the form 

1 - P(z7g Z7,) = 20^20^3 (1 - cj^cls) 

1/12 2 2 • 2/^^23-^\ fr,\ 

+ 4S12C13S13 sm ( ), (7) 

where the terms with Am^2 averaged as before. It can be justified because the argument 
of the sine term is of the order of 10 — 10^ with Am^ = 1 — lOOeV^, E= 4MeV, and L= 40m, 
the typical parameters of the Bugey experiment. The second term of (|^ may be neglected 
(as we will do) because sine-squared factor is ~ 10^^ for Am^ = 10~^eV^. 

We first examine the case that the atmospheric neutrino anomaly is attributed to the 
t'/x ~^ oscillations. Our discussion does not distinguish the types a and b until we address 
the constraint due to the double /5 decay. 

We demand, for consistency with the gross features of the LSND, the Bugey, and the 
atmospheric neutrino experiments, the following constraints: 

C?2C?3('5l2C23 + Ci2S23Sl3)^ =e< 10"^ (8) 



44(44 + 4) <^ = 2.5x10-2 



(9) 



(S12C23 + Ci2S23Si3)2(si2S23 - Ci2C23Sl3)^ < 0.1, (10) 
4C23S234(C12C23 " S12S23S13) (C12S23 + S12C23S13) ~ 1. (11) 

The constraints (|]) comes from the LSND experiment. We treat e as a small number of the 
order of ~ 10^^ or less. Our discussion will be insensitive to the number and we use it as a 
tentative guide when we address the consistency with other experiments. The equation (^ 
is due to the bound 1 — P{ve ^ '^e) < 5 % obtained in the Bugey experiment ||T^. It includes 



statistical and systematic uncertainties. The remaining two restrictions are from the features 
of the atmospheric neutrino data that the zenith-angle dependence is well described by an 
effective two-flavor-mixing ansatz with Am^ ~ lO-^eV^ and sin^ 2^ ~ 1. The constraint 
(|lO|) arises from a mild requirement that the first term of should be less than 0.2 so 
as not to disturb the effective two-flavor description. We emphasize that the constraints 
from the atmospheric neutrino data take the simple forms ([T0| ) and ( |TTD because of the mass 
hierarchy Am^2 ~ ^^13 ^ ^^23- 

We flrst notice that from the Bugey constraint (|^) that X = c\2c{^ must satisfy the 
inequality — X + 5 >Q. This inequality is so powerful that restricts the value of X into 
the two tiny regions < X < 5 and 1 — 5 < X < 1. On the other hand, we must have 
c^3 ~ 0(1) in order to satisfy the requirement (|n|). Thus, we have either c^2 — 5 or cfg — 1 
corresponding to the small-X and the large-X solutions, respectively. It is also required 
that C23S23 ~ I in order to maximize ([TT|). The small-X solution is then inconsistent with 
(ITTI). We end up with the unique solution 

(A) 4^4~e, 4-4-^ (12) 

where we have also utilized the LSND constraint (§) to push — 5 down to — e< 10^^. 
We have explicitly verifled that the allowed mixing pattern implied by (O) is physically 



unique throughout the varying mass hierarchies obtained by the cyclic permutations of 1-3 



of (I]), as it should be. Namely, the light "z/g" and the almost degenerate strongly mixed 
heavy "z/^" and "z/,-" for the type-a, and the heavy 'Ve" and the almost degenerate strongly 
mixed light "z/^" and "z/^" for the type-b cases. 

We now turn to the case that the atmospheric neutrino anomaly is caused by the z/^ — > z/g 
oscillation. In this case we replace the requirements (|T^) and ([IT]) by 

C?2C?3('5l2C23 + Ci2S23Sl3f < 0.1, (13) 

and 

- 4si2S23C?3Si3(ci2C23 - S12S23S13) ^ 1, (14) 

respectively. By the similar procedure one can show that the consistent solution of the 
requirements (H), (|^), (p!3|), and (|l^) is uniquely given by 

(B) 4^4^v^, c?3^s?3^1. (15) 

The solutions of the other type of mass hierarchies can be obtained by the similar manner 
and correspond to the redefinition of the mass eigenstates. The allowed mixing pattern is 
again physically unique: The light "z/,-" and the almost degenerate strongly mixed heavy 
"z/e" and "z/^" for the type-a, and the heavy "z/,-" and the almost degenerate strongly mixed 
light "z/e" and "z/^" for the type-b mass hierarchies. 

We note that the solutions (A) and (B) are subject to the additional constraints from 
other terrestrial experiments. While the solution (A) solves them automatically the nontriv- 
ial constraints arise for (B). In particular, the most stringent one comes from the Fermilab 
E531 experiment ||2^ for Am^ > 3eV^ and the z/^ disappearance experiment by the CDHS 
group ||2^ for Am^ < 3eV^. If we take the rate of the appearance events reported in at 
its face value the solution (B) may be excluded apart from tiny regions. To establish the 
rate, however, an additional run of the experiment as well as its careful analysis would be 
required. 

In the case of Majorana neutrinos a further constraint emerges from the non-observation 
of the neutrinoless double P decay. The quantity 
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< m^e >=^Vj\ Uej I '^mj (16) 



is constrained to be less than ~leV by the experiments |]T6[ where r]j = ±1 is the CP phase. 
Notice that we are working with the representation in which the mixing matrix is real under 
the assumption of CP invariance. 

Generally speaking, the constraint from the double f3 decay distinguishes between the 
type-a and the type-b mass hierarchies. In the type-a case there is a chance for cancellation 
between nearly degenerate two masses, but no chance in the type-b case because the heavy 
mass is carried by a unique mass eigenstate. 

New features, however, arise in our consistent solutions obtained above. We first discuss 
the case of atmospheric neutrino anomaly due to the — > z/^- oscillation. It can be shown 
that in the type-a mass pattern the double /? constraint is automatically satisfied because 
the heavy masses are always multipled by small angle factors. On the contrary, the angle 
factors in front of the unique heavy mass are always of the order of unity in the type-b mass 
hierarchy. Therefore, there is no consistent solution of the double /9-decay constraint for 
Majorana neutrinos in the type-b hierarchy. 

In the case of atmospheric neutrino anomaly due to the z/^ v^. oscillation, the situ- 
ation is somewhat different. In the type-a mass pattern there is a trouble because almost 



degenerate heavy masses are multiplied by 0(1) coefficients and a tuning, i.e., s^gC^s = 



13 



to better than 0.1, is required for cancellation in addition to the requirement of opposite 
CP parities. On the contrary, in the type-b hierarchy, there is no trouble with the double 
/3-decay constraint because the heavy mass is multiplied by small coefficients of the order of 

Thus, we have shown in this paper that the neutrino masses and the mixings are strongly 
constrained by the atmospheric and the terrestrial experiments under the assumption of 
DMMN suggested by the LSND experiments. The constraint is so severe that the mass 
and the mixing patterns are determined almost uniquely within the uncertainties of the 
neutrino types and the interpretations of the atmospheric neutrino anomaly. In the case of 



the Majorana neutrinos the additional constraint from the double (3 decay selects out the 
unique natural solution in each interpretation. 

Finally we give a few remarks. 
(1) Our analysis in this paper is less powerful in constraining the absolute values of the 
neutrino masses than in restricting the relative masses and the mixing angles. All the 
constraints would be cleared by, for example, the type-a solution with mi =6 eV, m2 =6.5 



eV, and = (6.5 + e) eV, which is also consistent with the direct measurements ||2T| and 



the cosmological considerations ||22 



(2) We have not performed a full three-flavor analysis of the atmospheric neutrino anomaly 
but relied on the effective two-flavor interpretation of either Uj. or channels. 

While the point deserves further study we are under a strong feeling that coexistence of the 
two channels does not spoil our solutions obtained in this paper. They certainly survive in 
the case of equal contributions of these two channels. 

This work is supported in part by Grant-in- Aid for Sciectific Research of the Ministry of 
Education, Science and Culture, #0560355. 

Note Added: After submitting the earlier version of this paper we became aware of two 
reports from the LSND group [^,^] with mutually conflicting conclusions. 
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